Loop quantization of spherically symmetric 
midi-superspaces : the interior problem 



Miguel Campiglia*, Rodolfo Gambini 1 ^ and Jorge Pullin** 

* Department of Physics The Pennsylvania State University, 
104 Davey Lab, University Park PA 16802 
^Institute de Fisica, Facultadde Ciencias, Igud 4225, esq. Mataojo, Montevideo, Uruguay. 
**2. Department of Physics and Astronomy, Louisiana State University, 
Baton Rouge, LA 70803-4001 

Abstract. We continue the study of spherically symmetric vacuum space-times in loop quantum 
gravity by treating the interior of a black hole. We start from a midi-superspace approach, but a 
simple gauge fixing leads to a Kantowski-Sachs form for the variables. We show that one can 
solve the quantum theory exactly in the (periodic) connection representation, including the inner 
product. The evolution can be solved exactly by de-parameterizing the theory and can be easily 
interpreted as a semi-classical evolution plus quantum corrections. A relational evolution can also 
be introduced in a precise manner, suggesting what may happen in situations where it is not possible 
to de-parameterize. We show that the singularity is replaced by a bounce at which quantum effects 
are important and that the extent of the region at the bounce where one departs from classical general 
relativity depends on the initial data. 
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INTRODUCTION 

In a previous paper [Q]] we have studied the quantization of spherically symmetric midi- 
superspaces in loop quantum gravity. We were able to treat the space-time exterior 
to a black hole in both the connection and the loop representation and recovered the 
quantization that Kucha? [2] had constructed using the traditional metric variables. The 
resulting quantum theories have wavefunctions that are functions of the mass of the 
space-time that do not evolve in time. The treatment being limited to the exterior, one 
could not probe the more interesting issue of what happens to the singularity in the 
quantum theory. In this paper we would like to address that issue. We will see that 
the same midi-superspace ansatz for the variables that we used for the exterior can 
accommodate the interior of a black hole. Through a simple gauge fixing one ends up 
with variables in Kantowski-Sachs form. We will construct the quantum theory in the 
connection representation. We will find that in the end one can actually solve the theory 
exactly, and it is easy to describe the evolution as a semi-classical one plus quantum 
corrections. We will see that generically if one starts with a state that approximates 
general relativity well at the cosmological horizon, quantum corrections eventually 
become important and a "bounce" replaces the cosmological singularity. The quantum 
solution therefore evolves past the point where one had the classical singularity and 
generically moves into another regime that approximates the classical theory well and 



develops a cosmological horizon. The Hamiltonian constraint that governs the evolution 
defines trajectories that depend on two constants of the motion. Even though in the 
classical theory one of them may be simply reinterpreted as a rescaling in the quantum 
regime both have physical consequences. Generically, the bounce is not symmetric: the 
extent of the regions where the quantum corrections are important is not the same before 
and after the bounce. However, the properties of the "inner" cosmological horizon are 
the same as the one one started with. If one were to naively match the resulting space- 
time to Schwarzschild one would have that the horizon after the bounce would have the 
same mass than the one of the horizon to the past of the singularity. We shall see that 
initial conditions can be set in such a way that the evolution is symmetric with respect to 
the region where the bounce occurs. The resulting symmetric wavefunctions will depend 
only of arbitrary functions of one parameter, the mass of the Black Hole. 



CLASSICAL SETUP 

We start with the same ansatz for the metric variables as in the exterior case. This 
is justified since the ansatz is good enough to encompass the interior Schwarzschild 
solution written in Kantowski-Sachs form. The kinematical setting for loop quantum 
gravity in spherically symmetric situations is well established and was discussed in 
detail by Bojowald and Swiderski [3] and can also be seen in our previous paper [Q]]. 
There is only one non-trivial spatial direction (the radial) which we call x since it is 
not necessarily parameterized by the usual radial coordinate. The canonical variables 
usual in loop quantum gravity are a set of triads Ef and SO(3) connections A' a ; after 
the imposition of spherical symmetry one is left with three pairs of canonical variables 
(t{,P^ jAqjEV ,A x ,E x ) (they are canonical up to factors involving the Immirzi parameter 
y). The variables and (p are angles transverse to the radial direction as in usual polar 
coordinates. Instead of using triads in the transverse directions, one introduces a "polar" 
set of variables £' <p , r\ and their canonical momenta. The relationship to more traditional 
metric variables is, 
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K XX = -sign(^)^^, Kee = -Vm^, (2) 

and the latter two are the components of the extrinsic curvature, y is the Immirzi 
parameter of loop quantum gravity. 

The Gauss law and diffeomorphism constraint read, 

G = P* 1 + (£*)' (3) 
D = />V +E ( f > A' (p -(E x )'A x . (4) 

and from now on we eliminate the rj variable and its conjugate momentum by solving 
Gauss' law and introduce K«, = A(p/2y and K x = (A x +7] r ) /2y. The diffeomorphism and 



Hamiltonian constraint now read, 
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This form of the Hamiltonian differs slightly from the one we used in yj] in the presence 
of the sign of the triad in the second term. In the exterior treatment of our previous paper 
the sign was constant and could be omitted, but it will change value when going from K x 
to — K x so we need to keep it here. We have also rescaled the lapse by a factor of sign of 
E x , so from now on N will be the ordinary lapse times the sign. The system has two pairs 
of canonical variables (K X ,E X ) and {K ( p,E <<> ) (we are taking Newton's constant G = 1). 
The constraints are first class. 

We now fix a gauge <j) = (E x )' = 0. Preservation of = in time partially determines 
one of the Lagrange multipliers, the lapse, 

= {(E x )',H} = (7) 
= (INKy)' = 0^N' = 0, (8) 

where we have used the diffeomorphism constraint, that together with <j) = implies 
that K'fp = 0. The latter is a secondary constraint that also has to be preserved in time. 
The preservation implies that (E^)' = 0, which is also a secondary constraint that is 
automatically preserved. Therefore all the variables are now independent of the radial 
coordinate x. The secondary constraints (E^)' = and KL = are second class and have 
to be imposed strongly. 

One is left with only one constraint, the Hamiltonian constraint, 



dxN + sg(E x )2K x K (p ^/\P\ + 
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In the above expression N.K^.E^ ,E X are only function of the evolution parameter, 
which we call X and are independent of the radial coordinate x. If one chooses K x in- 
dependent of the radial coordinate initially, one is left with a system where all canonical 
variables are functions of the evolution parameter, which we will call T, only. The inte- 
gral on x that appears in the constraint is over a finite interval Lq 0]. 

It is convenient to introduce variables that are more commonly used in loop quantum 
cosmology. The new variables are b,c,pb,p c which are canonical up to a factor of the 
Immirzi parameter and are defined as b = yKq,, c = jK x Lq p\, = E^Lq, p c = E x . With 
these variables, the Hamiltonian constraint has the form, 



p h 2bcp c lb 2 p b 
2 y 2 2 y 2 



where we have rescaled the constraint eliminating an overall factor 1 / ^fp~ c - Notice that 
a rescaling of the x coordinate also rescales Lq. This implies that a rescaling of x does 
not affect the value of pb and c. However, a change in Lq will rescale the value of pb and 
c. Therefore one can consider this dependence as an additional gauge freedom of the 
cosmological counterpart of the Schwarzschild solution. The evolution equations read, 
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c = — , (12) 

N 

Pb = —(2cpc + bp b ), (13) 
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2Nbpc 

Pc = —j—, (14) 

with N the rescaled lapse. 

The usual approach is to choose one of the variables as clock, for instance p c = £t (t) 2 
where £ = ±1 since an analysis of (|1 1H14|) shows that p c has a definite sign throughout 
the evolution. The classical evolution contains a singularity. Equation (fl4l) allows to 
determine the lapse, 

i(x)y 
bt{%) 

Substituting in (fTTI) one gets an ordinary differential equation for b 2 that can be imme- 
diately integrated, leading to, 
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where M is an integration constant. At the horizon one has that b = for the horizon to 
be isolated, so t = 2M there. From here we get, 



2M 

/> = eyj-l + ^, (17) 

where we have chosen the sign such that the original (un-rescaled) lapse remains positive 
when one tunnels through the singularity. 
Solving now (fT2l) we get, 

t(x) 2 



where Co is a constant. Using equation (1131) and the Hamiltonian constraint which 
implies that Pb = when b = 0, we get, 



p b = 2C y/t(T) (2M-t(T)). 



(19) 



The constant Co can be chosen by rescaling Lq, so we do so in such a way that 2Cq = 1 . 
If we choose the parameterization t = x the lapse is completely determined and takes 
the form, 

N= , 1 (20) 

2M _ i 

t 

This leaves us with the traditional line element for Kantowski-Sachs, 

ds 2 = - - l) l dt 2 + (If- - l) dx 2 + t 2 d& 2 . (21) 

We therefore see that we recover the traditional Kantowski-Sachs solution within 
the midisuperspace we started from. The system has only one constant of integration 
corresponding to one mechanical degree of freedom given by the value of M, just like in 
the exterior case. The slicing corresponding to the t = const, corresponds in the Kruskal 
diagram to hyperboloids r = const.. 



THE QUANTUM THEORY 
Quantization 

We can now proceed to "holonomize" the classical expressions found in order to 
carry out a loop quantization. The dynamical variables in the Hamiltonian are b,pb c 
and p c . Introducing a Bohr compactification for b and c one obtains a "loop quantum 
cosmology" which has been discussed by Ashtekar and Bojowald [4] and Modesto |@]. 
Instead of repeating their construction we will limit ourselves to introducing a technique 
that will be useful to discuss with more detail the issue of the bounce. More recently, 
Boehmer and Vandersloot |@] carried out a study of the holonomized semiclassical 
theory both with the \x and \x approaches. Their results for the \x case parallel the ones 
we discuss in the next subsection, although our approach is relational, whereas they 
work in a parameterized way. 

The Bohr compactification implies that the configuration variables take values in 
a compact space. The kinematical Hilbert space of wavefuctions is therefore given 
by periodic functions of the configuration variables. The "holonomized" version of 
the Hamiltonian constraint (fTOl) is an expression that is well defined as a quantum 
operator acting on such a Hilbert space. It therefore involves the configuration variables 
in periodic fashion. The expression depends on a parameter \i and reproduces the 
Hamiltonian (fTOl) in the limit ji — > 0. It takes the form, 

p b 2sin( i u% c sin( i uc) p h sm(^ib) 2 

2 n 2 f ^ 2 y 2 ' 

Upon quantization in the Bohr compactified space we keep ji finite, as is customary in 
loop quantum cosmology. The use of a fixed value of \l breaks the scale invariance of 
the Hamiltonian as it can be easily seen from (|22)) . 



In the Hilbert space considered the momentum variables are, 

p b w(b,c) = irf M ^^, 

p L V(b,c) = iyl^J^A, 



(23) 
(24) 



and on this space b and c are not well defined but periodic functions of them are. 

In order to quantize, it is better to rearrange the classical Hamiltonian constraint in a 
way that it is easy to "deparameterize", 
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which upon quantization yields the following form for the symmetric form of the "true" 
Hamiltonian, 

H Tm( W(b,c) = 



Z '^Planck7 
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(sin( M &) 2 + J uY)^^ (26) 
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Equations (1251261) can be solved using standard techniques for linear partial differ- 
ential equations that allow to reduce the system to a system of ordinary differential 
equations. The solution can be written in the form, 
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¥ (b, c)=A (b, c) exp ( -^S (b, c 
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with k a constant with dimensions of length squared, and where A is given by, 

A(b,c)=d, 



sin 2 (lib) +y 2 n 2 



(1 + 7 2 jU 2 ) \sm(/ib)\ 
where C\ is a normalization constant and S(b, c) is a solution of 
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(29) 



whose general solution is an arbitrary function S = &(w(b, c)) with w(b, c) a constant of 
the motion given by 



w(b,c) = In 



sin(/ic) 
cos(/i c) + 1 
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(30) 



Semiclassical approximation 



As in the usual eikonal approximation S contains the classical behavior of the solu- 
tions. Here, however, since we are working with the holonomized theory this will include 
the effective semiclassical dynamics of loop quantum gravity. This holonomized theory 
depends on an arbitrary parameter /i, as is usual in the minisuperspace context. If one 
wishes to identify the value of the parameter }i with similar concepts of the full theory, 
for instance, based on the quanta of area, the value of /i depends on Planck's constant. 
That is how one understands the classical holonomized theory as a semiclassical approx- 
imation to the dynamics of the model. The analysis can be done for any function 4>(w) 
and leads to solutions that differ by constant scalings, which can be reabsorbed in the 
constant k. We choose 4>(w) = w. Unlike the usual case in unconstrained systems, where 
it depends on two, here it depends on only one integration constant k. 

With the solution for S we can study the classical behavior of the system. In the usual 
Hamilton-Jacobi theory one starts by computing the canonical momenta of the original 
configuration variables, in this case b and c, 

4fcsin(/ifr)jU 

Ph ~ rV + l-cos^) 2 w 



k/i 



sin(jtic) 



(32) 



and the canonical momenta of the constants of integration, in this case only one, k, which 
are also constants of the motion, 

4tanh _1 ( cos(Mfc) 
p k = -ln( , Vv¥ ^ 7 (33) 
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We will now consider the initial conditions, which we will take at the cosmological 
horizon. We know that there one has b = (isolated horizon condition) and one identifies 
p c = 4M 2 so it can be isometric to the usual expressions in Schwarzschild. The initial 
value of c which we label cq is in principle arbitrary. Substituting in the equations (13 11321) 
we get Pb = initially and the constant k takes a value determined by cq and M and has 
dimensions of length squared. It is given by: 

AM 2 

k = sin(uco). (34) 

\x 

Finally, the equation for p^, which is independent of time, evaluated at different times, 
yields an equation from which we can obtain one of the configuration variables in terms 
of the other and the initial conditions. In particular, if we choose to solve for b, we get, 



b = —cos 



-l 
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(35) 



Substituting b into the expressions for pt, and p c we have all canonical variables as 
functions of c which operates as time variable. The analytic expressions are lengthy. To 
analyze their behavior it is best to resort to numerical evaluation. For this we need to 
choose numerical values for the various constants. For simplicity we take the Immirzi 
parameter to be unity 7=1, which is close to the value that stems from black hole 
entropy. The parameter jl is arbitrary for the classical theory, although usual loop 
quantum cosmology arguments lead to /i = y/3/4. However this leads to a "bounce" 
happening far from the Planck radius. This led in loop quantum cosmology to consider 
a variable value of }l. We will not pursue that route here. We will just consider /i as 
a free (but small) parameter. In the numerical computations shown here we have taken 
jl = 0.02. 

The problem has two constants of the motion. One of them is the initial value of c = cq. 
The other is the value of the "mass" — if one is viewing the problem as a Schwarzschild 
interior — , which is determined by the value of k. When the problem is being considered 
as a Schwarzschild interior, it is natural to rescale cq by varying Lq as we noted before. 
In the holonomized theory, however, this rescaling freedom is not present at the quantum 
level due to the finite value of the "area quantum" /i, and different values of cq lead to 
different behaviors, in particular, at the bounce. We shall analyze this effect and show 
that a preferred value is naturally selected once /i is fixed. 

An interesting check is to note that in the limit fl — > the semi-classical solution of 
the holonomized quantum theory recovers the traditional Kantowski-Sachs space-time 
of general relativity, as expected. In particular one has that the independent constant Co 
of Kantowski-Sachs is related with the initial constant cq as Co = AcqM/ y. 

We have evaluated numerically the exact semi-classical solutions given before in order 
to study their properties. In figure \T\ we show the volume as a function of the time 
variable c. One sees a bounce occurring where the singularity would have appeared 
in general relativity. The volume does not vanish and the co-triad does not diverge. 
The volume goes to zero at the Kantowski-Sachs horizon, which occurs at the past an 
future horizons where b = 0. We have chosen b = at the initial data, where we chose 
c = cq = 10~ 6 (the value cq = yields an ill-defined evolution). The volume then also 
vanishes at the final point, at c = nj\i — cq, irrespective of the initial value of cq. 

In figure [2] we show the volume as a function of the "radial" variable r 2 = p c . The 
volume is multi-valued and it corresponds to the contracting and expanding phases. The 
rate at which the volume contracts depends on the value of cq. We have chosen it here in 
such a way that the contracting and expanding rates are very close to each other. This, 
however, is a choice. In particular, the region where the holonomized theory departs 
from general relativity will change depending on such choice. If one chooses, as we 
did, a "symmetric" evolution, the region of departure is confined near the bounce in 
a symmetric way. Other choices can lead to the region extending far away from the 
bounce either into the future or the past of it. It is remarkable that no matter what choice 
of co initially, the system will eventually go past the bounce and return to a regime 
where the holonomized theory approximates general relativity and a future horizon will 
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FIGURE 1. The volume as a function of the time variable c. We see a bounce occurring in the place 
where usual Kantowski-Sachs would have a singularity (left panel). The right panel shows the volume as 
a function of time close to the Kantowski-Sachs horizon (b = 0), where the volume vanishes. Calculations 
use M = 5 x 10 5 M Planck . 



develop with the same mass value. That is, the presence of both horizons and the bounce 
is robust. How the bounce occurs and when the holonomized theory starts and stops 
its departure from general relativity depends on the initial data. This is a reflection of 
the instabilities that were noticed in the recursion relations of the loop representation 
treatment of these models 
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FIGURE 2. The volume as a function of the "radial" time variable r 2 = p c . It can be seen that it 
contracts till it achieves a minimum value, it then expands at the same expansion rate as it contracted. 
This is achieved if one chooses the initial data for cq at the horizon in such a way that the evolution 
is symmetric at the "bounce". This, in particular implies that the region where the holonomized theory 
departs from general relativity is minimized. For other choices, the region of departure is larger. 




Figure [3] shows the sine of fib versus the time variable c. The sine of fib divided by 
fl is the variable in the holonomized theory that replaces the connection Aa, in the usual 
theory. Both theories agree when the sine is small. We see that indeed the sine is small 
close to the horizons and maximal at the bounce. The time variable c chosen is such that 
the region close to the "bounce" measured in terms of the radial variable is achieved 
very fast and lasts for a long "time", therefore the sine of \lb is large through most of 
the evolution. One can see that when measured in terms of the volume the connection 
departs from it classical value close to bounce. 
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FIGURE 3. The sine of fib vs c. When the sine is close to one the holonomized theory departs from 
general relativity significantly. We see that both theories agree at the initial and final moments, when 
the cosmological horizons are present. The diagram is symmetric in spite of its appearance, it is hard 
for the plotting program to numerically get to the points near the final horizon, since the points gather 
exponentially. 



Quantum corrections 

We have discussed the classical approximation. Let us now proceed to consider the 
quantum corrections, which involve the prefactor A(b,c) in equation (1271) . We start 
by noting that the proposed solution for the wavefunction satisfies the boundary con- 
ditions implied by the kinematical Hilbert space. That is, v I / (— 7t/fl,c) = *P(7r/jU,c), 
^(b, —n/fl) — *F(&,7r/]Li). We have required that the boundary conditions be periodic 
in c and b with period 2ltJ fl. It turns out that the only region of phase space of impor- 
tance is between and it/ fl. In the other three regions one obtains the same |*F| 2 and 
the physical predictions are the same. It might appear surprising to impose boundary 
conditions in the variable c, which we are using as "time". In principle one would think 
that specifying "future" boundary values would violate causality. Remarkably, this is not 
the case. The value of cq is arbitrary and the final result of the evolution is in no way 
determined by the periodicity imposed. 




The dynamical Hilbert space is therefore given by arbitrary functions of the two 
constants of the motion cq,M. As the Schwarzschild external solution depends only of 
M one can restrict the Hilbert space to the symmetric semi-classical solution around the 
bounce that are the ones whose behavior departs from the classical only close to radial 
distances of the order of Planck scale. Given /i this is equivalent to choose cq such that 
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and A = y 1 + y 2 /! 2 . Going to the full quantum theory, however, will not allow to choose 
precisely the values of c$,M that correspond to the symmetric bounce, and one will have 
a spread implying that the region where the bounce occurs will be larger than the one 
we consider here in the symmetric case. 

Finally, one can show that the evolution in c is unitary. In fact one can see that Hj me 
is self-adjoint by using the following technique [9]: Check the dimensionality of the 
two subspaces: Jtf + = ker(Hj me + i) and = ker(//xme _ They are self-adjoint if 
and only if both spaces have the same dimensionality. In this case the wavefunctions *F 
multiplied by the exponential factor 



exp I =1=4 tanlf 



_j / cos {lib) sin (/xc) 
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(38) 



respectively belong to each kernel and are normalizable in the kinematical space. An 
important observation is in order here, the elements of the physical space ^(b.c) are 
also normalizable with the inner product of the kinematical space and therefore they 
can be also used to define a relational evolution in terms of conditional probabilities, 
as for instance in the proposal of Page and Wootters [1(1 U]. A detailed study of these 
definition of the evolution that do not require a deparametrization of the constraint will 
be given elsewhere. 

We have therefore carried out a quantization of the interior of the Schwarzschild 
space-time using loop quantum gravity techniques. We have shown how the singularity 
is replaced by a bounce, and what conditions are needed for the bounce to occur in a 
regime of Planck energy. We have also outlined how one would construct the quantum 
theory for the model. It is interesting to compare the interior and exterior treatments. In 
the latter (and also in the complete space-time treatment of Kuchaf [2]) one is left with an 
quantum theory in which wavefunctions are arbitrary functions mass of the space-time, 
which is conserved. Here, while treating the interior as a cosmology, we are left with 
arbitrary functions of an observable, which evaluated on the initial data is determined 
entirely by cq. This is the variable conjugate to p c , whose initial value is associated 



with the mass of the Schwarzschild space-time. Therefore the two pictures are clearly 
reconciled. 
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